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■ Abstract 

r~| ! The beam transfer matrix, often called the ABCD matrix, is a 

two- by- two matrix with unit determinant, and with three independent 
parameters. It is noted that this matrix cannot always be diagonal- 
ized. It can however be brought by rotation to a matrix with equal 
diagonal elements. This equi-diagonal matrix can then be squeeze- 
^ ' transformed to a rotation, to a squeeze, or to one of the two shear 

. matrices. It is noted that these one-parameter matrices constitute the 

basic elements of the Wigner's little group for space-time symmetries 
^SJ ' of elementary particles. Thus every ABCD matrix can be written 

If^ . as a similarity transformation of one of the Wigner matrices, while 

I the transformation matrix is a rotation preceded by a squeeze. This 

mathematical property enables us to compute scattering processes in 
periodic systems. Laser cavities and multilayer optics are discussed 
in detail. For both cases, it is shown possible to write the one-cycle 
^ ' transfer matrix as a similarity transformation of one of the Wigner 

■ matrices. It is thus possible to calculate the ABCD matrix for an 
arbitrary number of cycles. 
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1 Introduction 



Without too much exaggeration, modem physics consists of harmonic oscilla- 
tors and two-by-two matrices. Ray optics is mostly the physics of two-by-two 
matrices. The beam transfer matrix, usually called the ABCD matrix, is a 
two-by-two matrix with unit determinant and with three independent pa- 
rameters. 

In mathematics, the group of matrices with this property is called the 
two-dimensional symplectic group or 5*^(2). Indeed, the 5*^(2) group has 
been playing a pivotal role both in classical and quantum optics. It is the 
basic language for lens optics [21 E]; multilayer optics ^E], laser cavities [0], 
as well as squeezed states of light in quantum optics |71 |B] ■ 

This group is locally isomorphic to the group of Lorentz transformations 
applicable to one time dimension and two space dimensions. This Lorentz 
group consists of three-by-three matrices. Since the algebra of two-by-two 
matrices is much simpler than those of three-by-three, the Sp{2) group pro- 
vides a very convenient computational tool for Lorentz transformations in 
physics [9]. 

While our understanding of the Pauli matrices and the rotation group is 
thorough, we are not yet completely familiar with the matrices of the Sp{2) 
group pm [TT] . Do we know how to diagonalize the ABCD matrix? We 
shall point out in this paper that, unlike the case of the rotation matrices, 
it is not always possible to digonalize the ABCD matrix by a similarity 
transformation. We shall show however that every ABCD matrix can be 
rotated to a matrix with equal diagonal elements. 

We shall then show that this equi-diagonal matrix can be squeezed to 
one of the four one-parameter matrices consisting rotation, squeeze, and two 
shear matrices. They take the form 



cos 6* — sin 6*\ / cosh A* sinh A* 
sin 9* cos 6* J ' \ sinh A* cosh A* 

1 -7*\ fl 

1 y' V7* 1 
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(1) 



One of the two triangular matrices in the above expression can be obtained 
from the other by rotation, but it is more convenient to use both of them. 

It should be noted that these four matrices form the four basic matrices of 
Wigner's little group for space-time symmetries for elementary particles 
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. We shall therefore call them Wigner matrices, and use W collectively 
for these four matrices. 

These matrices have equal diagonal elements, but not all of them are 
diagonalizable. The first matrix is a rotation matrix, whose properties are 
well known. However, its eigenvalues are complex, and therefore it cannot 
be diagonalized if all the matrices in the group are to be real. The second 
matrix is symmetric and can be diagonalized by a rotation. The third and 
fourth matrices are triangular and cannot be diagonalized. 

On the other hand, these matrices share one convenient feature with the 
diagonal matrix. Namely, for the first matrix, 

cos 9i — sin Oi\ / cos 62 — sin 62 
sin 61 cos Oi J \ sin 62 cos 62 



cos (^1 + 62) - sin (^1 + 62) 
sin {Oi + 62) cos (^1 + 62) 



(2) 



(3) 



•cos^ -sin^\^ _ /cos(A^^) -sin(A^^)\ 
sin^ cose ) ~ Uin(iV^) cos(iV^) j' 

The multiplication of two matrices results in the addition of parameters. This 
was called "slide-rule" property in reference [5] , but it is more appropriate to 
call it "logarithmic property." The remaining three matrices in equation ([1]) 
have this logarithmic property. 

We shall show in this paper that the ABCD matrix can be brought to one 
of the Wigner matrices by a similarity transformation, and that the similarity 
transformation is a rotation followed by a squeeze. 

This mathematical result can be useful in calculating the ABCD matrix 
for periodic systems. In this paper, we study laser cavities and multilayer 
optics in detail. In both cases, we reduce the multi-cycles system into one 
cycle. Both subjects have been extensively discussed in the literature. 

Recently, the present authors used the method of Wigner's little group to 
calculate the ABCD matrix for laser cavities [2]. However, their cycle had 
to start from the midway between the two mirrors. In this paper, our cycle 
will start at an arbitrary point. 

In 2003, Georgieva and Kim used a method based on the Lorentz group to 
study multilayer optics [5] . These authors also end up with the inconvenience 
of starting their cycle from the midway in one of the layers. We eliminate 
this inconvenience. 
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In Sec. m the eigenvalues of the ABCD matrix is discussed in detail. They 
can be real, complex, and one. Not all of them can be brought to a diagonal 
form, especially if we restrict the matrices to be real. It is noted however that 
every ABCD matrix can be brought to a form with equal diagonal elements 
by a rotation. In Sec. [3|, we construct a similarity transformation which will 
bring the ABCD matrix into one of the four Wigner matrices. It is shown 
that this matrix is a rotation matrix followed by a squeeze matrix. It is also 
shown that these four different Wigner matrices can be combined into one 
analytical expression. 

The present formalism is applicable to one-dimensional periodic systems. 
In Sec. m we study laser cavities in detail. In Sec. we show how useful 
Wigner matrices are in computing the scattering matrix for multilayer optics. 



2 Eigenvalues of the ABCD matrix 

Let us start with the two- by-two matrix 

A B 
C D 



(4) 



where the elements A,B,C,D are real. The determinant of this matrix is 
one and thus {AD - BC) = 1. 

The eigenvalue equation for this two-by-two matrix is only a quadratic 
equation, and it is always soluble. The eigenvalue equation becomes 

- {A + D)E + AD - BC = 0, (5) 

and the eigenvalues are 

1 



E± = -i^A + D±^{A- Df + 45Cj . (6) 

If the quantity inside the square root sign is negative and {A+D)/2 is smaller 
than one, then the eigenvalues become 

E± = exp(±ir), (7) 

with 

cosr = ^^. (8) 
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If the quantity inside the square root sign is positive, and {A+D)/2 is greater 
than one. The eigenvalues are 



E± = exp (±A*), 



(9) 




(10) 



1. 



(11) 



These are all well known, and the purpose of this section is not to repeat 
these trivial procedures. 

Although we are used to diagonal matrices with eigenvalues as diagonal 
elements, it is very cumbersome to deal with the above three separate cases. 
In order to find a common ground, we propose to use the matrix of the form 
with equal diagonal elements: 



and propose to bring the ABCD matrix of equation (jl]) to this form by 
a similarity transformation. Since the trace is preserved under similarity 
t r ansf ormat ions , 



Since the determinant of the ABCD matrix is one. 



In this paper, we started with the ABCD matrix with three independent 
parameters. Because of the condition on the diagonal elements, the JFCJ 
matrix has two independent parameters. If the off diagonal elements F and 
C have opposite signs, it can be written as 




(12) 



-FC= 1. 



(14) 




(15) 



If F and G have the same sign, it can take the form 



( 



cosh A* 
e~'^ sinh A* 



e'' sinh A* 
cos A* 



(16) 
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If F or G vanishes, we can write the JFGJ matrix as 

,0 1 j' U-V J- ^^^^ 

Let us go back to equation ([1]), the above matrices are similarity trans- 
formations of the Wigner matrices with the transformation matrix 



e-^'^ 



(18) 



which expands one axis and contracts the other in the two-dimensional space. 
This is a squeeze transformation. In general, symmetric two-by-two matrices 
perform squeeze transformations. Indeed, the matrices in equations (|T5l) to 
f|T7|) are "squeezed" Wigner matrices. 



3 Similarity Transformations 

In order to find a similarity transformation which will bring the ABCD 
matrix of equation (j4]) to the equi-diagonal matrix of equation (fT2|l . let us 
try a rotation 

cos(5/2) sin(5/2)\M 5Wcos(5/2) - sin(5/2) \ 
-sin((5/2) cohi5/2))\C D)\sin{5/2) cos((5/2); 

We can calculate the angle 5 first. The result is 

tan5 = ^^. (20) 

Using this angle, we can then compute F and G. They are 
_ 1 



(19) 



B-C+J{A-Df + {B + CY 



G=]^{c-B + ^{A~DY + {B + Cyy (21) 

If J is smaller than one, we can write J = cos 6*, and use equation f[T^ 
for the JFGJ matrix, with 

cos9 = — - — , 



^ c-B-^jA-py + jB + cy 

C-B + J(A-Dy + {B + Cf 
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(23) 



and the Wigner matrix should take the form 

cos 9* — sin 6* 
sin 9* cos 6* 

which is one of the four matrices given in equation ([1]). 
If J is greater than 1, we should use equation (fT6|l . with 

^ + D 
cosh A = — ^ — , 

, B-C + J {A - Df + (5 + Cf ^ ^ 

C-B + MA-DY + {B + Cy 



The Wigner matrix in this case is 



/ cosh A* sin A* 
V cosh A* sinh A* 



(25) 



which is also given in equation ([T]). 

If J = 1, the diagonal elements are one, and the off-diagonal element F 
or G has to vanish. As a consequence, the ABCD matrix becomes 

I ^7^) {c-B ?)' ^26) 

and 

C-5 = eV, or C-B = e-'^-i\ (27) 

which are similarity transformations of the triangularr matrices given in equa- 
tion ^. 

Therefore, ABCD is a similarity transformation of one of the Wigner 
matrices given in equation ([1]). We thus write 

ABCD = SWS-\ (28) 

where W is the Wigner matrix given in equation ([T]), and the similarity 
transformation matrix is 

/cos(5/2) -sin(5/2)We^/2 \ 

I sin(5/2) cos(5/2) A e'^/^ J ■ K^^) 

This is a rotation matrix preceded by a squeeze matrix. 
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If there are four different Wigner matrices, calculations become very cum- 
bersome in practical applications. There are no problems if only one Wigner 
matrix, such as that of equation fl23|) . is applicable throughout the problem 
as in the case of laser cavities discussed in Sec. HI 

If two or more forms are involved in one problem, as in the case of multi- 
layer optics discussed in Sec. \5\ the problem could be uncontrollable. Thus, 
we are led to consider one mathematical expression for all four Wigner ma- 
trices. For this purpose, let us consider the form 

cos(^/2) -sin(^/2)\ /coshA sinh A \ / cos(^/2) -sm{9/2)\ , . 
sm{9/2) cos(^/2) j UinhA cosh A j V sin(^/2) cos(0/2) j ' ^ ^ 

This is a special case of the Bargmann decomposition [13], and its physical 
interpretation in special relativity is discussed in the literature In this 
paper, we consider this form purely for mathematical convenience for dis- 
cussing periodic systems. After multiplication, the above expression can be 
compressed into 

/ (cosh A) cos ^ — (cosh A) sin 6' + sinh A' 
V (cosh A) sin 6' + sinh A (cosh A) cos ^ 



(31) 



with equal diagonal elements. This matrix has two independent parameters. 
Indeed, this is another form of the JFGJ matrix with 

J = (cosh A) cos6', 

F = — (cosh A) sin 9 + sinh A, 

a = (cosh A) sin 9 + sinh A. (32) 



When A = 0, it becomes 



while it takes the form 



cos 9 — sin 9 
sin 9 cos 9 



(33) 
(34) 



/ cosh A sinh A 

V sinh A cosh A , 

when ^ = 0. 

When tanh A = ± sin 9, the equi-diagonal matrix of equation (1311) becomes 

1 0\ /I -2 sinh 

,2 sinh 7] ij' VO 1 

respectively. 



(35) 
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4 Laser Cavities 



The laser cavity consists of two concave mirrors separated by distance s. The 
mirror matrix takes the form 



1 0' 
-2/R 1, 

where R is the radius of the concave mirror. The separation matrix is 

1 s 
1 



(36) 



(37) 



If we start the cycle from one of the two mirrors one complete cycle consists 

-2/R l)(o l)(-2% l)(o l)- ^^^^ 
If we start the at the position d from the mirror, then the complete cycle 
becomes 

1 d\ f 1 0\ /I s\ / 1 0\ /I s-d\ 

(40) 



lJ\-2/R lj\0 l)\-2/R lAO 1 



smce 



1 s\ _ (1 s-d\ fl d' 
,0 ij ~ VO 1 j U 1 
Thus, one complete cycle consists of two repeated applications of the half 
cycle, which can be written as 

which becomes 

(l-2d/R {l-2d/R){s-d)+d\ , . 

^ ~ [ -2/R 1-2(3- d)/R ) ■ ^^'^> 

One complete cycle then becomes L^. There is however inconvenient fea- 
ture of this expression is that the off-diagonal elements have different dimen- 
sions, while the diagonal elements are dimensionless. In order to deal with 
this problem, we write this expression as a similarity transformation 

\/l-2d/R {l-2d/R){s- d)/s + d/s\ (1/^ 
l/^/sJ\ -2s/R l-2{s-d)/R )\ y/s 

(43) 
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The elements of the middle matrix are now dimensionless. Let us now use 
the notations a = d/ s and b = 2/R. Then we can work with the normalized 
L, which takes the form 



1 — ab 1 — ab{l — a) 
-b 1-6(1 -a) 



(44) 



where s = 1, according to the normalization defined in equation (H3|l . We 
can now bring this form into a equi-diagonal matrix. We can write this 
expression in the form of equation (fT9|) with 



tan 5 



b-2ab 



1 — b — ab(l — a) ' 
F = i (^(1 + 6) - a6(l - a) + ^ (b - 2aby + [1 - b {1 + a - a^)f^ , 

G = i (^ab{l -a)-{l + b) + ^{b- 2abf + [1 - 6 (1 + a - a^)f^ . (45) 

We can now write the ABCD matrix for one complete cycle as a similarity 
transformation 

-- SR{Ae*)S-^, 



and 



with 



SR{2e*)S-^ 

R(2e*) -- 



COS 9* — sin^^* 
sin 9* cos 9* 



cosr = 1 - -. 

2 



(46) 
(47) 

(48) 



The similarity transformation matrix 5* is in equation fl29l) and takes the 
form 

S -- 



e''/2 cos(5/2) sin(5/2) \ 

6^/2 sin(5/2) 6-^/2 cos(5/2) J ' 



with 



(49) 
(50) 



where F, G, and 5 are given in equation fH51) . 

Thus, for cycles, the ABCD matrix becomes 



2N 



SR{29*)S-^ = SR{4:N9*)S-\ 



(51) 
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where we have used the logarithmic property of the Wigner matrix. 

When the cycle starts from the midpoint in the cavity, a = 1/2, and the 
half-cycle matrix becomes 



1-6/2 1-6/4 
-6 1-6/2 



(52) 



(53) 



and 6 becomes zero. We do not need the rotation matrix, and this matrix 
can be written as 

cos^* e''sin^* 
-e"''sin6'* cos6'* 

with 

/)* 1 ^ 2« 4-6 

cost^ =1 , e ' = — — . 

2' 46 

The matrix L can of course be written as a similarity transformation 



3^/2 \ / cos^* sinrWe-''/2 g 
e-'i/y [-sine* cos 6* J \ e'i/\ 

For the laser consisting of cycles, the ABCD matrix becomes 



(54) 



T2N_(e'^'^ \(cos{Ne*) sin(iVr) We-''/2 q \ , . 
~\Q -e-^l^)\-sm{Ne*) cos {N9*) J [ -e^/^' ^^^^ 

2N^f cos(iVn e^sin(iVr)\ 

\,-e~'^sm{N9*) cos (NO*) 7' ^ ^ 

This is the result we obtained in our previous paper on laser cavities [B]. 
The point of this paper is we can start the cycle from an arbitrary point, by 
introducing the parameter 6. 

When the cycle begins from one of the lenses, a = 0, and the L matrix 
becomes 

1 1 
-6 1-6. 



(57) 



This expression can first be brought to the equi-diagonal form by a rotation 
matrix with the rotation angle 

tan5 = ( — ^ ) , (58) 
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according to equation (I2UI) . The equi-diagonal matrix then takes the form of 
equation (1551) or equation with 



cose* = 1 - -. 

2' 



, (1 + 6) + J62 + (1_5)2 

e^^ = - - V ^ (59) 

(1 + 6) - + (1 - 6)2 



Indeed the L matrix of equation ( l57l) is a similarity transformation of the 
Wigner matrix 

cos 6* sin 6** 
— sin 9* cos 61* 



(60) 



with the transformation matrix of the form given in equation (129|) . which is 
a rotation preceded by a squeeze. Then it is straight-forward to calculate the 
ABCD matrix for cycles. 



5 Multilayer Optics 

In multilayer optics, we are led to consider two beams moving in opposite 
directions pj. One is the incident beam and the other is the reflected beam. 
We can represent them as a two component column matrix 



E+e 



ikx 
-ikx 



(61) 



where the upper and lower components correspond to the incoming and re- 
flected beams respectively. For a given frequency, the wave number depends 
on the index of the refraction. Thus, if the beam travels along the distance 
d, the column matrix should be multiplied by the two-by-two matrix [1] 



Q 

e-*"/2 



(62) 



where a/2 = kd. Thus, the propagation matrices for two different media can 
be represented by 



g*m/2 Q X /gm2/2 g 

Q g-iai/2 ) > ( Q g-»a2/2 



(63) 
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respectively, with ai/2 = kid and 02/2 = k2d. 

If the beam propagates along the first medium and meets the boundary 
at the second medium, it will be partially reflected and partially transmitted. 
The boundary matrix is [H H] 

cosh(/i/2) sinh(/i/2)\ , . 

sinh(/i/2) cosh(/x/2) / ' ^ ' 

with 

cosh(;u/2) = l/ti2, sinh(/i/2) = ri2/ti2, (65) 

where ti2 and ri2 are the transmission and reflection coefficients respectively, 
and they satisfy (rjg + ^12) ~ 1- The boundary matrix for the second to first 
medium is the inverse of the above matrix and can be written as 



(66) 



cosh(/i/2) - sinh(^/2) \ 
-sinh(/i/2) cosh(/i/2) J 

Thus one complete cycle, starting from the second medium, consists of 

cosh(/i/2) sinh(/i/2) \ /e*"i/2 g 
sinh(/i/2) cosh(^/2) j V e-^^^/^ 

/ cosh(/i/2) -sinh(/x/2)\ n 

^-sinh(^/2) cosh(/i/2) A g-^^^/a ) • l^g 



This expression can be unitarily transformed to [16] 

\ /cos(ai/2) - sin(ai/2) \ /g-^/^ 
sin(ai/2) cos(a;i/2) /V e^/^ 

/cos(a2/2) -sin(a2/2)\ , s 

Uin(a2/2) cos(a2/2) J' ^ > 

The first three matrices in the above expression can be compressed to 

cos(ai/2) -e'^ sin(ai/2) \ 
e-^sin(ai/2) cos(ai/2) j' 

and thus, according to equation fl3T|) and equation fl33|) . to 

cos(6'i) cosh A — sin(6'i) cosh A + sinh A 
sin(6'i) cosh A + sinh A cos(6'i) cosh A 



(69) 



(70) 
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with 

cosh A = (cosh /i) \Jl — cos^ (cn /2) tanh^ /x, 

cos(q;i/2) ,^ . 

cos0i = 1 (71) 

(cosh/i) A/l — cos^(ai/2) tanh /i 



Then the matrix M of equation ( 1681) can be written as 

/cos(0i/2) -sin(0i/2)\ /coshA sinhA 
Uin(^i/2) cos(0i/2) j UinhA cosh A 

/cos(e2/2) -sin(e2/2)\ , . 

Uin(e2/2) cos(^2/2) )' ^ ' 



with 



02 = Oi + a2. 



We can now write this expression as a similarity transformation of a 
matrix with equal diagonal elements. Explicitly, 

cos(5/2) - sin((5/2) \ / J F\ / cos(5/2) sin(5/2) \ 
sin(5/2) cos{5/2))[g j)[~sm{5/2) cos(5/2) j ' ^ ^ 



where 



with 



J F\ _ /cos(^/2) -sin(e/2)\ /coshA sinh A 
G J J ~ \sm{e/2) cos(^/2) j UinhA cosh A 

cos(^/2) -sin(e/2)\ 

sin(^/2) cos(^/2) r ^ ^ 



^'^-02), 9 = 1 (9, + 92) . (75) 



2 ' " 2 

The above three matrices can be compressed to one equi-diagonal matrix 



cosh A cos 9 — cosh A sin 9 + sinh A 
cosh A sin 9 + sinh A cosh A cos 9 



(76) 



This matrix can now be written as a similarity transformation of one of the 
four Wigner matrices. 
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When the off-diagonal elements have opposite signs, we can write this as 



e°"/^ \ f cos LU —smuj\ ( e '^/^ 
e 



"/"^ Vsincx; coscx; M e'"/^ 



(77) 



with 



cosh oj = cosh A cos 9, 

2^ cosh A sin 6^ — sinh A ^^^^ 
cosh A sin 6 + sinh A 



When the off-diagonal elements have the same sign, we should write 

(79) 



with 



J F\ _ /e'"/2 \ /coshx sinh x \ / e'''/^ q 
G e""^/ V gjni^ ^ coshxJ V e"/^ J ' 



cosh X = cosh A cos 9, 



sinh A — cosh A sin 6* , , 

e = . (80) 

cosh A sin 6 + sinh A 

If one of the off-diagonal elements vanishes, the matrix takes the form 

/I -2sinhA\ / 1 0\ 

U I ) UsinhA ij- 

Now the diagonalization of the above expression is straight-forward ac- 
cording to the mathematical tool we have developed in this paper. So is the 
computation of 

M^. (82) 

The computation of this process was started in reference [5]. There, the cycle 
had to start from the midway in one of the media, but no explanation was 
given why. In this paper, by introducing the R{S) matrix, we can start from 
an arbitrary point. Indeed, this is a physical interpretation of this rotation 
matrix. 

When the cycle starts from an arbitrary point in the medium, we can 
adjust the value of the rotation angle 6 as we did in Sec. HI 
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Concluding Remarks 



In dealing with a two-by-two matrix, we are accustomed to think it can be 
diagonahzed and can be brought to the diagonal matrix by a rotation. This 
paper shows this assumption is not always true, because the ABCD matrix 
is not always a rotation matrix. 

It has been shown that the ABCD matrix can be brought to one of 
the four one-parameter Wigner matrices through a similarity transforma- 
tion, and the similarity transformation is a rotation followed by a squeeze 
transformation. 

These one-parameter Wigner matrices have the logarithmic property which 
allows us to calculate repeated applications by multiplying the parameter by 
an integer. This property is transmitted to its similarity transformation, 
since 



This mathematical technique is applied to laser cavities and multilayer 
systems. However, it is clear from the literature that the mathematical result 
stated in this paper is a result of our efforts to understand physical systems. 

This mathematical instrument is applicable to other periodic systems in 
physics, such as one-dimensional scattering problems in quantum mechan- 
ics [HllTB], especially in condensed-matter physics. We can also broaden our 
scope to look into applications in space-time symmetries of elementary parti- 
cles in view of the fact that the Wigner matrices used in this paper are from 
Wigner's 1939 paper on symmetries in the Lorentz-covariant world [T2| [19]. 

For periodic systems, many authors used different approaches. Sanchez- 
Soto and his co-authors used the conformal representation of the Sp{2) group 
to attack the problem (HED]- This is possible due to the fact that the trans- 
formations of this group can be translated into conformal transformations as 
noted by Bargmann [14j. 

Another interesting approach to the periodic system is to use mathemat- 
ical induction [21]. It is possible to assume first that the ABCD matrix is 
known for cycles, and then compute the system for + 1 cycles. 




(83) 



16 



References 

[1] Azzam RAM and Bashara I 1977 Ellipsometry and Polarized Light 
(North-Holland, Amsterdam) 

[2] Ba§kal S and Kim Y S 2007 Phys. Rev. E 63 056606-056611 

[3] Ba§kal S and Kim Y S Phys. Rev. E 67 056601-056608 

[4] Monzon J J and Sanchez-Soto L L 2000 J. Opt. Soc. Am. A 17 1475-1481 

[5] Georgieva E and Kim Y S 2003 Phys. Rev. E 68 026606-026612 

[6] Ba§kal S and Kim Y S 2002 Phys. Rev. E 66 026604-026609 

[7] Yuen H P 1976 Phys. Rev. A 13 2226-2243 

[8] Kim Y S and Noz M E 1991 Phase Space Picture of Quantum Mechanics 
(World Scientific, Singapore) Chap 5 

[9] Kim Y S and Noz M E 1983 Am. J. Phys. 51 368-375 

[10] Gillemin V and Sternberg S 1984 Symplectic techniques in physics (Cam- 
bridge University Press, Cambridge) Chap 1. 

[11] Lang S 1985 SL2{R) (Springer- Verlag, Heidelberg) Chap 1 and 5 

[12] Wigner E 1939 Ann. Math. 40 149-204 

[13] Kim Y S and Noz M E 1986 Theory and Applications of the Poincare 
Group (Reidel, Dordrecht) Chaps 3, 4 and 9. 

[14] Bargmann V 1947 Ann. Math. 48 568-641 

[15] Han D and Kim Y S 1988 Phys. Rev. A 37 4494-4496 

[16] Georgieva E and Kim Y S 2001 Phys. Rev. E 64 026602-026607 

[17] Sprung D W L, Wu Hu, and Martorcll J 1993 Am. J. Phys. 61 1118-1124 

[18] Grifiiths D J and Steinke C A 2001 Am. J. Phys. 69 137-154 

[19] Han D, Kim Y S, and Son D 1986 J. Math. Phys. 27 2228-2235 



17 



[20] Sanchez-Soto L L, Carinena J F, Barriuso A G, and Monzon L L 2005 
Eur. J. Phys. 26 469-480, and the references contained in this paper. 

[21] Kildemo M, Hunderi O, and Drevillon B, 1997 J. Opt. Soc. Am. A 14 
931-939 



18 



